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Three AGN Close To The Effective Eddington Limit For Dusty Gas 



CO 

o 

CN 
X> 

IX, 

oo 
CN 

X 

Q-i! 

6 

CO 

o3 



> 
ON 

o 
o 
q 

cn 
o 
cn 



R.V. Vasudevan 1 *, A.C. Fabian 2 , R. F. Mushotzky 1 , M. Melendez 1 , L. M. Winter 3 
and M. L. Trippe 1 

1 Department of Astronomy, University of Maryland, College Park MD, 20742, USA 

2 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, 

3 Atmospheric and Environmental Research, 131 Hartwell Avenue, Lexington, MA, USA 

4 March 2013 



ABSTRACT 

The Effective Eddington Limit for dusty gas surrounding AGN is lower than the canonical 
Eddington limit for hydrogen gas. Previous results from the Swift/BAT 9-month catalogue 
suggested that in the overwhelming majority of local AGN, the dusty absorbing gas is below 
this Effective Eddington limit, implying that radiation pressure is insufficient to blow away 
the absorbing clouds. We present an analysis of three objects from that sample which were 
found to be close to the Effective Eddington limit (NGC454, 2MASX J03565655-4041453 
and XSS J05054-2348), using newly obtained XMM-Newton data. We use the X-ray data 
to better constrain the absorbing column density, and supplement them with XMM optical 
monitor (OM) data, infrared Spitzer and Herschel data where available to construct a broad- 
band spectral energy distribution to estimate refined bolometric luminosities and Eddington 
ratios for these three objects. The new XMM-Newton observations show all three objects 
moving away from the region expected for short-lived absorption in the Ah — AEdd plane into 
the 'long-lived absorption' region. We find our conclusions robust to different methods for 
estimating the bolometric luminosity and Eddington ratio. Interestingly, 2MASX J03565655- 
4041453 and XSS J05054-2348 now exhibit complex X-ray spectra, at variance with previous 
analyses of their Swift/XRT data. We find evidence for absorption variability in NGC 454 and 
2MASX J03565655-4041453, perhaps implying that although the radiation pressure from the 
central engine is insufficient to cause clearly detectable outflows, it may cause absorption 
variations over longer timescales. However, more robust black hole mass estimates would 
improve the accuracy of the Eddington ratio estimates for these objects. 

Key words: black hole physics - galaxies: active - quasars: general - galaxies: Seyfert 



1 INTRODUCTION 

The tight relation between the mass of the central black hole of 
a galaxy and the mass or velocity dispersion of the surrounding 
stellar bulge points to a coupling or feedback between the black 
hole and the host galaxy. The accretion power produced by the 
growth of the black hole exceeds the binding energy of the galaxy 
by one to two orders of magnitude so some form of feedback is 
expected. What is unclear however is how it works, and direct ev- 
idence is elusive and difficult to obtain. Radiation pressure from 
the Active Galactic Nucleus (AGN) will have little effect on the 
stars of the surrounding bulge and only act on the gas. Catching 
and recogni sing an object whe re the gas is being expelled has not 
been easy. Fabian et all J2009h find clear evidence for the effect of 
radiation pressure from the central AGN on surrounding dusty gas, 
using the most complete hard X-ray selected sample o f AGN avail- 
able at the time, the Swift/BAT 9-month catalogue iTueller et alj 
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2008). This catalogue consists of 97 low redshift AGN with mea- 
sured column densities 7V~h and inferred black h ole masses. They 
use the X-ray spectral fits of IWinter etail (§009) (W09 hereafter) 
to measure the amount of gas, parameterised as the column den- 
sity A*h, and use X-ray luminosities in conjunction with a bolo- 
metric correction (to provide the full ionising continuum luminos- 
ity) and a black hole mass to parameterise the effect of radiation 
pressure, as the Eddington ratio AEdd- The Eddington ratio AEdd 
is the ratio of the bolometric luminosity of an AGN (Lbd) t0 me 
Eddington luminosity, or 'Eddington Limit' for its black hole mass, 
L Edd = 1.3 x 10 38 (A/ BH /M o ). 

The canonical Eddington limit for Hydrogen and Helium re- 
lies only on electron scattering. When metals and dust are mixed in 
cold gas, the cross section to radiation becomes much larger, due 
to the high absorption cross section of dust in the UV, so the effec- 
tive Eddington limit decreases. The deviation from the canonical 
Eddington limit depends on the spectral energy distribution of the 
central AGN and the column density of the absorbing clouds. Fig.Q] 
shows the variation of the effective Eddington limit as a function of 
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Nn dFabian et al. 2008), The exact locus of this effective Eddington 
limit also depends on details such as the dust-to-gas ratio of the sur- 
rounding material, but at high column densities this has little effect. 
The radius of the clouds from the black hole (BH) also alters their 
effective Eddington limit since the stars inward from the cloud will 
also exert inward gravitational force, but even if the stars contribute 
as much mass as the central BH, this will only increase the Edding- 
ton limit by a fac tor of two. This is re presented by the blue dashed 
line in Fig.0 In lFabian etal] b009l) . results from the SwiftfB AT 
9-month catalogue (wit h Eddington ratios ca lculated using bolo- 
metric corrections from IVasudevan & Fabian 20071) show a clear 
avoidance of the 'forbidden' region in the Nn — AEdd plane, point- 
ing to the pronounced influence of the AGN on its environment. 
Sources well below the Eddington limit do not have sufficient ra- 
diation pressure to exert significant influence on their host galaxies 
and so their absorption is stable. Absorbed sources at low column 
densities may appear above their effective Eddington limit because 
of the presence of dust lanes in the galaxy (i.e. the mass interior to 
the position of the clouds is much greater than the central black hole 
mass). The sources on the boundary, near or above the effective Ed- 
dington limit, are of particular interest since they must be exhibit- 
ing an episode of pronounced interaction between the AGN and 
the absorbing gas. The one source which lies well into the forbid- 
den region is MCG-05-23-16, whi ch is known to hav e a complex 
and variable warm absorber (e.g. iReeves et al.ll2007h and a pos- 
sible variable, high-ionisation, high- velocity outflow dBraito et alj 
120071) . thus implying it is near the Eddington limit. For warm ab- 
sorbers, the effect radiation pressure is reduced. Detailed further 
study of other near-effective Eddington limited/super-effective Ed- 
dington sources such as these could be extremely illuminating in 
understanding the role of AGN feedback in shaping the host galaxy. 

There are three ne ar-effective Eddington limited sources from 
the iFabian etak I J2009h study which previously only had poorer- 
quality Swift/XKT data - NGC 454, 2MASX J03565655-4041453 
and XSS J05054-2348 - the XRT data were insufficient to con- 
strain the detailed nature of the absorption in these sources. We 
successfully proposed to obtain XMM-Newton observations for 
these three luminous absorbed AGN. These sources have complex 
optical spectra; for NGC 454, inspection of optical spectra from 
the REOSC Spectrograph at C ASLEO 2.15 m Ritchey-C hretien 
telescope, San Juan, Argentina (Donz elli & Pasto riza 2000) shows 
a dual system with a Seyfert 2 (eastern source, NGC 454E) 
and a starburst galaxy (western source, NGC 454W); 6dF spec- 
tra for 2MASX J03565655-4041453 and XSS J05054-2348 re- 
veal that these two objects can be classified as a starburst galaxy 
and a Seyfert 2 respectively (u sing the line ratio diagnostics of 
IVeilleux & Osterbrocld[T987l and lKewlev et alj|2006t) . 

These three objects were previously observed for total periods 
of 20 ks (NGC 454, 345 counts), 12.6 ks (2MASX J03565655- 
4041453, 790 counts) and 3.2 ks (XSS J05054-2348, 389 counts) 
using XRT and the analyses of their spectra are presented in W09. 
For NGC 454, the XRT spectrum can be fit by a partial covering or 
double power law model giving a strong initial indication that com- 
plex absorption is at work, and the other two objects are fit by sim- 
ple power-laws. However, given the low signal-to-noise ratio of the 
Swift/XRT data, such models may not be physically appropriate to 
fully describe these systems and the acquisition of far superior qual- 
ity XMM spectra allows us to fit more detailed models. The aim 
with this study is to constrain whether outflows are present and con- 
strain the nature of the absorption (cold or warm/ionized). One of 
these XMM datasets ( for NGC 454) has been analysed in detail by 
iMarchese et alj ( |2012|) . and we build upon their conclusions in light 
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Figure 1. The Effective Eddington Limit paradigm: types of absorption ex- 
pected if AGN are located in different areas of the Nn — ^Edd plane. The 
Effective Eddington Limit, calculated using CLOUDY simulations, is shown 
by the blue solid curve (assuming standard interstellar medium grain prop- 
erties), and differs from the canonical Eddington Limit (Agdd = !)■ The 
dashed blue curve shows a factor of 2 increase in the effective Eddington 
limit due to the mass of intervening stars equal to the black hole mass itself 
(therefore providing a stronger gravitational counter-balance to radiation 
pressure for dusty gas clouds sufficiently far out from the nucleus). AGN 
below this limit are expected to have long-lived absorption. AGN above this 
limit, in the 'forbidden region' are expected to have short-lived absorption 
due to significant radiation pressure on the absorbing clouds, and should 
therefore not stay in this region for prolonged periods. Column densities 
log(Nn) < 21.5 may be due to dust lanes in the host galaxy. 



of the Effective Eddington Limit paradigm here. The XMM spectra 
allow more accurate determination of their X-ray luminosities and 
column densities, refining their position on the Ah — AEdd plane. 
Determinations of source variability between different observation 
epochs will be valuable in constraining how much these sources 
move around in this plane. For these objects, the high column den- 
sity makes it difficult to recover the bolometric luminosity from 
optic al-to-X-ray SED integration, so we employ archival Spitzer 
and Herschel data where available to provide an alternative esti- 
mate, for calculation of Eddington ratios. Additionally, we require 
good estimates of the black hole mass to determine the Edding- 
ton ratio. We employ published black hole mass estimates using 
th e black hole mass-K-ban d host galaxy bulge luminosity relation 
of iMarconi & Hun 3 d2003l) fr om th e literature (although we note 
that the later study ofjGraham 2007 provides an updated, more ro- 
bust version of this relation for future determinations). The approx- 
imate method outlined in W09 uses the archival 2MASS extended 
source magnitudes as an estimate of the total galaxy light (includ- 
ing the nucleus) and subtracts the published point source magni- 
tudes _Jor_each_source toremove the nuclear co ntribution. How- 
ever, Vasudevan et al. (2009) (V09 hereafter) and I Vasudevan et alj 
(2010) (V 10 hereafter) find that this over-estimates the bulge lumi- 
nosity considerably, since the bulges of most of these galaxies are 
unresolved at the redshifts typical for the catalogue. Therefore, they 
use the fainter point source magnitudes as the total 'bulge plus nu- 
cleus' light, and use infrared SED templates to predict the level 
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of nuclear contribution that needs to be removed in the K-band 
to get the bulge contribution alone. We use the refined method of 
V09/V10 for our black hole mass estimates whenever possible. It is 
not possible to obtain virial black h ole mass estimates u sing optical 
broad lines for these sources (e.g jBaskin & Laorl l2005) since the 
optical spectra for these sources indicate either Seyfert 2 or star- 
burst types and do not display broad lines. 



2 DATA SOURCES AND REDUCTION 

The details of the XMM data used are presented in Table Q] We 
reduce the raw data using the XMM-Science Analysis Software 
(SAS) using the EPCHAIN and EMCHAIN tools as described in 
IVasudevan & Fabianl d2009l) . including correcting for flaring by in- 
specting when the background lightcurve exhibited flaring episodes 
that were a significant fraction (up to a fifth) of the source flux 
and excluding those parts of the observation from further process- 
ing. We extract the PN, MOS1 and MOS2 spectra using 36 arcsec 
source regions (corresponding to an encircled energy fraction of 
~ 85 per cent for PN and MOS instruments) and multiple back- 
ground regions near the source, using the standard procedures de- 
tailed in the XMM-SAS manual. We group the final X-ray spectra 
with a minimum of 20 counts per bin, to allow \ 2 statistics to be 
used in fitting the data in XSPEC. We extract the optical-UV OM 
data using the OMICHAIN pipeline, checking for a good match be- 
tween the optical source and the X-ray source position in each filter, 
and de-redden the OM fluxes for Galactic du st extinction us i ng on- 
line dust extinct ion maps generated f rom the lSchlegel et alj dl 
dataH using the lCardellietai1dl989ri Galactic extinction curve. 

We supplement the optical, UV and X-ray data with infrared 
data from Spitzer and Herschel when available. The Spitzer/IKS 
spectra are availab le in the Cornell Atlas of Spitzer/IRS sources 
dLebouteiller et all 1201 lh . and we find data for NGC 454 and 
XSS J05054-2348 in these archives. We only find Herschel data 
for XSS J05054-234 8; this objec t was mapped wit h the PACS 
dPoglitsch etai1l20ld) and SPIRE iGriffin et al]|20ld) instrument s 
on board the Herschel Space Observatory dPilbratt et al.l l2010l) . 
Data were taken simultaneously with the PACS blue and red chan- 
nels, centered on 70 and 160um. SPIRE observations were taken si- 
multaneously at 250, 350, and 500um. The PACS and SPIRE data 
were reduced in standard way to level 1 using HIPE 9.0.0 dOtq 
2010) and the final maps were made using the SCANAMORPHOS 
software (Roussel 2012, submitted). A detailed description of the 
Herschel data reduction process will be presented in Melendez et 
al. (in prep). 



3 X-RAY SPECTRAL FITS (XMM PN/MOS AND 
SWIFT/BAT) 

We fit the combined PN, MOS 1 and MOS2 spectra along with the lat- 
est 58-month Swift/BAT spectra, testing a suite of heuristic models 
using the XSPEC software to identify the best-fitting mode l, using 
the systematic method outlined in IVasudevan et all d2013h (Table 
2). In addition to attempting model fits using the model combina- 
tions in that paper (which include tests for Iron lines at 6.4 keV, 
soft excesses below 1 keV and edges due to highly ionised oxy- 
gen transitions at 0.73 and 0.87 keV), we also attempt other fits 
including a double power-law models (a direct component at high 
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Figure 2. NGC 454E: best fit model: an absorbed (Galactic and intrinsic 
absorption) double-power law with an Iron-line at 6.4 keV, reduced chi- 
squared: 1.23. The upper panel shows the model fit to the data (in counts 
space) with a ratio plot, and the lower panel shows the unfolded spectrum 
in v F„ units. 



energies and scattered component at lower energies with the same 
photon index T, analogous to partially-covering absorption) and in- 
vestigate whether low-energy residuals can be modelled as due to 
hot, diffuse gas using the MEKAL model. 

We carefully order the different model fits based on their re- 
duced chi-squared goodness-of-fit (x 2 /d.o.f). The best-fit model 
results for all three objects are presented in Table |2] where the er- 
rors on all quoted parameters are 90 per cent confidence intervals 
provided by the 'error' command in XSPEC. 



3.1 NGC 454 

In the XMM (PN and MOS) images, we see one prominent point 
source surrounded by a few faint sources, which can be ascribed to 
a num ber of different optical counterparts. The work of Llohanssonl 
dl988l) identifies optical components to NGC 454 identified as east- 
ern (E), western (W), and three southern components (SE, S and 
SW). The prominent X-ray source could potentially encompass 
emission from both the E and W components (the two merging 
galaxies) , but analy s is of the line ratios for the two sources pre- 
sented in Johansson ( 1988) would strongly suggest that the eastern 
component is responsible for the AGN emission, and is identified 
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AGN Redshift RA Dec Observation Observation Observation Total counts Usable % of obs. 

ID date duration (ks) (PN and MOS) (without flaring) 

NGC454 0.012158 18.605 -55.397 0605090301 2009-11-05 29.9 4548 98 

2MASX J03565655-4041453 0.074782 59.235 -40.696 0605090201 2009-09-12 39.5 9614 84 

XSS J05054-2348 0.035041 76.441 -23.854 0605090101 2009-08-06 30.9 25241 98 



Table 1. Table of XMM-Newton observations used for each object. 



AGN 


T(l) 


JVjf* (2) 
(10 22 cm- 2 ) 


E FeK (3) 

(keV) 


EKcK (4) 
(eV) 


T unabs fc\ 
^2-10keV *- J > 
(lO^ergs -1 ) 


Scattered 
fraction (%) (6) 


X 2 /d.o.f. 


NGC 454E 

2MASX J03565655-4041453 
XSS J05054-2348 


1.87 ±0.06 

1 63+ 013 
J ~ DO -0.06 
i ci +0.05 
± ' ol -0.04 


7 1+0.6 

9.9 ±0.3 


6.34 ±0.05 
6.43 ± 0.06 
6.37 ±0.03 


1»3_ 102 


0.1624 ±0.0003 
3.54 ±0.02 
3.031 ± 0.003 


0.6 
1.0 
0.8 


1.23 (384.2/313) 
0.98 (681.0/694) 
1.06(1286.0/1215) 



Table 2. Fit results for complex power-law best-fits to XMM-Newton and Swift/BAT data. (1): photon index. (2): absorbing column density. (3), (4): En- 
ergy/equivalent width of iron line as modelled by a redshifted gaussian (ZGAUSS) component. (6) scattered fraction: the ratio of normalisations of power-law 
components (low-energy scattered power law normalisation / high energy direct power law normalisation). This may represent an upper limit on the true 
scattered fraciton, since the soft emission may be contaminated by thermal, non-AGN emission. 



in ljohanssonl Jl988h as a red elliptical galaxy in comparison to the 
western component, a blue, irregular (possibly disk) galaxy. Analy- 
sis of Chandra images for NGC 454 (Mushotzky and Koss, private 
communication) also shows that the hard emission (above 2 keV) 
is centred on the eastern source, confirming its AGN nature. We 
henceforth refer to this object as NGC 454E for clarity, to empha- 
sise that this is the source responsible for the AGN emission of 
interest here. 

The best-fit model is an absorbed double power-law 
with coupled photon indices (partial covering absorption) with 
an Iron line at 6.4 keV. In XSPEC, the model prescrip- 
tion is TBabs*zTBabs(powerlaw + zTBabs(powerlaw + 
ZGAUSS)), where the TBABS component accounts for fixed Galac- 
tic absorption (determined using the NH utility from FLTOOLS) and 
the ZTBABS components account for absorption intrinsic to the 
source; the TBABS model is adopted since it is the most up-to-date 
absorption model currently available in XSPEC. The fit parameters 
are detailed in Table [2] We also attempt reflection and ionised ab- 
sorption models, but do not find any significant improvement in the 
fit. There are some prominent residuals at low energies (~1 keV) 
which could be identified as emission from hot, diffuse gas indica- 
tive of star formation contamination to the AGN spectrum from the 
western component. The fit is marginally improved by inclusion of 
a MEKAL model to model these low-energy residuals (an improve- 
ment of A\ 2 ~ 7). There are also unmodelled residuals near the 
Iron K-alpha line ( ~ 5 — 7 keV) which ha ve been attributed to an 
ionised absorber in iMarchese et al.l d2012l) (see also %. It . but we 
do not attempt such a detailed fit here. 

W09 find a column of 15.9 x 10 23 cm" 2 from XRT spectra, 
in comparison to 36.8+ 2 '3 x 10 22 cm -2 found here using XMM. 
The variation of a factor of ~ 2 suggest this is a 'changing-look' 
AGN (e.g. see Risaliti et al .120021) , a scenario discussed extensively 
in lMarchese et al.l d2012l) by comparing the XMM-Newton dataset 
used here with Suzaku and Swift/XRT data. The intrinsic luminos- 
ity found here is very similar (within 13 per cent) of the value found 
in W09. The IMarchese et alj d2012l) analysis finds that the 2-10 
keV intrinsic luminosity varies by a factor of ~3-4 between XMM 
and Suzaku observations along with a very pronounced variation 



The BAT spectrum for this source is consistent with a power- 
law, and does not display any of the curvature expected from Comp- 
ton reflection. To check this, we fit the combined XMM and BAT 
data with a simple TBABS(PCFABS(PEXRAV)) model. We ignore 
the Iron line region 5.5-7.5 keV, which is not modelled by PEXRAV. 
This model includes Galactic absorption, intrinsic (partially cov- 
ered) absorption and a reflection continuum parameterised by the 
photon index T, the reflection amplitude R and the fold energy 
£fold (in keV). We freeze the other parameters in the PEXRAV 
model at their default values (abundances are solar, and the cosine 
of the inclination angle is fixed at 0.45), and check the validity of 
the fit by plotting contours of R against Sfoid • We allow the rela- 
tive normalisations of the BAT and PN to vary with respect to each 
other, since the BAT spectra are not simultaneous with the XMM 
spectra. This fit yields a reflection amplitude consistent with R = 
(R < 1.25), a fold energy £Jfoid = 56l2o° keV and a photon index 



of r = 1.29 



+0.09 



, much harder than the index of 1.87 obtained 



in absorption (between 



and ~ 2 x 10 /4 cm- 



using their modelling approach outlined in Table 1 of their paper). 



using the double power-law fit. In these fits, there can be degen- 
eracies when tryin g to independently con strain R, F and -Efoid all 
independently (see lVasudevan et alj|2013h : for example, the reflec- 
tion and photon index can be correlated over a wide range of values 
if the high-energy cut-off is not well constr ained. Additionall y, val- 
ues of T < 1.5 are not typical for AGN dCorral et aLll201ll) and 
are energetically difficult to obtain in s tandard Compton-scat tering 
paradigms for the coronal emission (Zdziarski et al.l [T990t) . We 
check the robustness of our reflection analysis by freezing the pho- 
ton index at the value 1.87 in the PEXRAV fit and find that this 
results in a weakly constrained fold energy (with only a lower limit 
of Efoid > 106 keV obtained from the error analysis) with a reflec- 
tion fraction 13.9i"f 'f , indicating that the uncertain relative normal- 
isation of the BAT and XMM data makes it difficult to constrain 
reflection in this source. The large value of reflection here would 
indicate a complex geometry; R = 1 typically indicates a covering 
solid angle for the reflecting material of 2n (representing, for ex- 
ample, an isotropic source above a disc), but values larger than this 
may be indicative of general relativistic effects such as light bend- 
ing, or complex absorption, neither o f which can be constrained 
with these data. IMarchese et al.l d2012) also attempt a reflection fit 
to Suzaku and Swift-BAT data, but separate out the direct and re- 
flected components rather than using a single PEXRAV model as 
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Figure 3. 2MASX J03565655-4041453: Best fit modekan absorbed (Galac- 
tic and intrinsic absorption) double-power law with an Iron-line at 6.4 keV, 
reduced chi-squared 0.98. The upper panel shows the model fit to the data 
with a ratio plot, and the lower panel shows the unfolded spectrum in v F„ 
units. 

done here. They find a reflection fraction of R ~ 1 (i.e. a reflector 
covering 2n steradians), consistent with what we find here, and no- 
tably show a very good agreement between the Suzaku Hard X-ray 
Detector (HXD) data and the Swift-BAT data (Fig. 4 of their paper). 

3.2 2MASX J03565655-4041453 

The best-fitting model to the combined XMM and BAT 
data is again a double power-law with an Iron line at 
6.4 keV (TBabs*zTBabs(powerlaw + zTBabs(powerlaw 
+ ZGAUSS))), as for NGC 454E (see Table|2]for details). The ab- 
sorbing column found here is 7.ligl4 x 10 22 cm~ 2 which again 
represents an increase in a factor of ~ 2 over the absorbing column 
found in W09 using an XRT spectrum. The intrinsic luminosity 
found here is about half that found in W09, suggesting significant 
intrinsic variability as well as a substantial variation in the absorb- 
ing column density. Additionally, the XMM data reveal a complex 
spectrum (resembling partially-covered absorption) as opposed to 
the simple absorbed power-law shape found using the XRT data in 
W09. 

We also perform a reflection analysis here using the XSPEC 
model combination TBABS(ZPCFABS(PEXRAV)), as done for NGC 



XSS J05054-2348 
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Figure 4. XSS J05054-2348: Best fit model: an absorbed (Galactic and in- 
trinsic absorption) double-power law with an Iron-line at 6.4 keV, reduced 
chi-squared: 1 .06. The upper panel shows the model fit to the data with a 
ratio plot, and the lower panel shows the unfolded spectrum in v F„ units. 

454E. The shape of the BAT spectrum could be suggestive of reflec- 
tion, and the best fit reflection parameters are well constrained: the 
reflection fraction is R — 3.4.1?, g, with a weakly constrained fold 
energy above 200 keV, and a photon index T = 1.84±g;?|, consis- 
tent (within lrr errors) with the value 1.63 obtained from the dou- 
ble power-law fit. If we explicitly freeze the photon index at 1.63, 
we find a reflection fraction consistent with zero R < 1.32 and 
a weakly constrained fold energy outside the BAT band (_Bf id > 
182 keV). Although the reflection fraction is consistent with zero 
(within errors) if we freeze the photon index, the better constrained 
photon index for this source may suggest that reflection is present 
on some level, potentially with a high reflection fraction. A high ef- 
fective reflection fraction could indicate light bending or extremely 
complex absorption in this source. 

3.3 XSS J05054-2348 

A double power-law with an Iron line at 6.4 keV 
(TBabs*zTBabs(powerlaw + zTBabs(powerlaw + 
ZGAUSS))) is again the best-fitting model for the combined XMM 
and BAT data (see Table O for details. We note the presence of 
some structure in the data below 2 keV, which can be well fit by a 
MEKAL model, improving the overall fit significantly ( A\ 2 ~ 30), 
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but do not discuss this further here. The absorbing column found 
here is 9.9 ± 0.3 x 10 22 cm -2 , representing an increase in a 
factor of ~ 1.5 over the absorbing column found in W09 using 
an XRT spectrum. The intrinsic luminosity found here is about 75 
per cent of that found in W09, suggesting both moderate intrinsic 
luminosity and column density variation. Here also, the XMM 
data reveal a complex spectrum (resembling partially-covered 
absorption) as opposed to the simple absorbed power-law shape 
found using the XRT data in W09. 

A reflection analysis, using the same approach as for the 
other objects, yields the best fit reflection parameters R — 



(R < 0.05), with a fold energy E {o 



44_g , and a photon 



index T = 0.95±g;g|, much harder than the value 1.51 obtained 
from the double power-law fit. This extreme photon index could 
indicate very complex absorption in this source. Freezing the pho- 
ton index at 1.51 obtained from the double power-law fit yields 
R = 11.9±o.g with EfoLd = 157^52 keV, so there is considerable 
systematic uncertainty in the reflection parameters depending on 
the assumptions adopted. 



3.4 Comparison with archival XRT and Suzaku data 

These three sources have all been studied in W09 using Swift/XKT 
data. We perform a brief comparison of our new XMM data with 
the XRT data used in W09, and fit the best fit model obtained 
for XMM+BAT data to XRT+BAT data. For all three objects, the 
absorbed double power-law fit yields very similar parameters to 
the W09 fits to the XRT data alone (within errors), although the 
XRT data do not have sufficient counts to constrain the Iron line 
parameters. NGC 454E has too few counts below 1 keV in the 
XRT spectrum to obtain a well-constrained photon index, and fit- 
ting the XRT data alone, W09 find V = 2.81± u ;jjf, whereas we 
find T = 1.91±o when including the BAT data. We show a 
comparison between the new XMM datasets and the XRT datasets 
in Fig. [5] from fitting the double power-law model jointly to both 
datasets (we do this to purely show the degree of variation between 
the spectra when fit with a common model, and do not inteipret the 
model fit results since it is not robust to draw physical conclusions 
from jointly fitting observations taken at different epochs). 

An archival S uzaku observa t ion of XSS J05054-2348 has been 
studied in detail bv lEguchi et al.l ( 120091) . They fit a complex, multi- 
absorber model to the combined Suzaku XIS, PIN/HXD and Swift- 
BAT data to obtain the absorption level in different components. It 
is not straightforward to perform a direct comparison with their 
work due to the different modelling approach, but one of their 
absorbing components has a similar column density to the value 
found here (their value is 5.87 ± 0.08 x 10 22 cm -2 , compared to 
9.9 ± 0.3 x 10 22 cm~ 2 here, although the hard component in their 
study has a higher absoiption of 29.3 ± 6.2 x 10 22 cm -2 ). Regard- 
less of the model fit, the shape of the Suzaku spectrum is qualita- 
tively very similar to that of the XMM spectrum taken three years 
later (see Fig. [6); this is intriguing given that the Suzaku light curve 
from 2 006 shows very lit tle variability compared to other sources 
in the lEguchi et alj fc009l) study. 



4 EDDINGTON RATIO ESTIMATES 

In order to refine the position of these three objects on the Nn — 
^Edd plot, we also require m ore accurate determ inations of the Ed- 
dington ratios. Previously, in lFabianetalJ ( l2009l) , Eddington ratios 
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Figure 5. Joint fits to both XMM (PN) and archival XRT data for the three 
objects. The higher line below 10 keV in each plot is the XMM data (binned 
for 2MASX J03565655-4041453 due to high counts); the lower line below 
10 keV is the XRT data (used in W09) and above 10 keV, we show the BAT 
data. 



were calculated using a bolometric correction applied to the X- 
ray luminosity (to estimate the bolometric luminosity Lboi). and 
div iding by the Eddin gton luminosity estimates calculated using 
the I Winter e t al. (2009) black hole masses. Here we refine this ap- 
proach by constructing spectral energy distributions (SEDs) from 
infrared to X-rays, and consider the best approach for determining 
Lbol in these three AGN. We make use of the XMM optical mon- 
itor (OM) data to construct the optical-to-Xray SED, but find that 
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XSS J05054-2348 




Figure 6. Suzaku, XMM and BAT data for XSS J05054-2348, jointly fit 
with a double power-law with absorption and an Iron line. The top line 
below 10 keV is the XMM data, the lower line below 10 keV is the Suzaku 
XISO data, the higher line above 10 keV is the Suzaku HXD data, and the 
lower line above 10 keV is the BAT data. 



all three of our objects exhibit very reddened optical-UV spectra, 
probably due to significant dust in our line of sight or host-galaxy 
contamination. This makes it difficult to recover the intrinsic ion- 
ising continuum responsible for exerting radiation pressure on the 
surrounding absorption. If the significant X-ray absorption by neu- 
tral gas found in fj3] is coupled with dust as expected in AGN, 
this would also suggest that the red optical-UV SED shapes are 
likely due to significant dust reddening. We therefore use two dif- 
ferent techniques to estimate the bolometric luminosity: 1) we at- 
tempt to de-redden the optical-UV continuum for dust by includ- 
ing intrinsic reddening corrections in the optical-UV regime and fit 
an accretion disc and power-law model combination to to recover 
the full intrinsic lumi nosity from optical-to-X-ray wavelengths 
( Vasu devan et~ai1l2009h . and 2) we use the reprocessed infrared ra- 
diation (from mid-to-far infrared wavelengths, 8-1000 /jm) to es- 
timate the portion of the intrinsic luminosity that is absorbed by 
circumnuclear dusty gas and re-emitted in the infrared and add it 
to the 'remainder' of the intrinsic emission that appears at high en- 
ergies above ~1 keV dLusso et"aT]|2012l IVasudevan et alj|2010l 
IPozzi et alj|2007h . We compare and contrast the results from both 
approaches, and determine whether the refined Eddington ratios for 
these objects still place t hem close to t he Eff ective Eddington limit 
for dusty gas as found in Fabian et al. (20091). 



4.1 Bolometric luminosity 

For using the optical-UV SED for determining Lboi, we use the 
Optical Monitor (OM) data to construct SEDs in conjunction with 
the PN data for the thr ee objects, and determin e bolometric lumi- 
nosities as detailed in [Vasudevan et alj J2009h . using black hole 
mass estimates to constrain the accretion disc model fit in the 
optical-UV. We ext ract the magnitud e s from the OM point sources 
as det ailed in e.g. Brocksopp et al. (2006) and I Vasudevan et al] 
( 2009) and construct the broad-band SED along with the PN data. 
We also include the Swift/BAT data in the 14-195 keV energy 
range. For the model fit, we use the best-fit found for the X-ray data 
alone as a starting point, and add in a DISKPN component (a mul- 
ticolour black-body accretion disc model) with intrinsic reddening 
incorporated using a ZDUST component to model the optical-UV 
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Figure 7. NGC 454E broad-band SED. Large black filled points represent 
OM and PN data. Above 10 — 3 keV, thin solid lines represent the model fits 
with absorption, heavy solid lines are absorption-corrected fits. The long- 
dashed-short-dashed line linking 13.6 e V to 1 keV is the slop e predicted 
using mid-IR high-ionisation line ratios dMelendez et alj|201lL see [|Hfor 
details on this method). Below 10~ 3 keV, small points represent Spitzer 
data, the dotted line represents the nuclear IR emission component and the 
solid line represents the combined nuclear and host-galaxy emission (in this 
case dominated by the nuclear emission alone). 
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Figure 8. 2MASX J03565655-4041453 broad-band SED. Black filled 
points: OM and PN data. Thin dotted lines represent the model fits with 
absorption, heavy solid lines are absorption-corrected fits. The dot-dashed 
lines linking 13.6 eV to 1 keV show the range of 13.6 eV emission that can 
be obtained using the slopes predicted (1 .5 < a < 2.0 for F v oc v~ a ) 
using mid-IR high-ionisation line ratios jMelendez et alj|201lL see ^5]for 
details on this method) found for the BAT catalogue. For clarity, we only 
plot the model fits that assume the W09 mass, since the fit using the V10 
mass is almost identical. 
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Figure 9. XSS J05054-2348 broad-band SED. Large black filled points rep- 
resent OM and PN data. Above 10 — 3 keV, thin solid lines represent the 
model fits with absorption, heavy solid lines are absorption-corrected fits. 
The dot-dashed line linking 13.6 eV to 1 keV is the minimal slope 1 . 5 pre- 
dicted using mid-IR high-ionisation line ratios ( Melen dez et al .2011, see 
for details on this method) found for the BAT catalogue. Below 10 — 3 keV, 
small points represent Spitzer data, diamonds represent Herschel data, the 
dotted line represents the nuclear IR emission component, the dashed line 
represents the IR host-galaxy component and the solid line represents the 
combined nuclear and host-galaxy emission. For clarity, we only plot the 
fits that assume the W09 mass, since the fit using the V10 mass is almost 
identical. 



accretion disc emission along with intrinsic reddening (assuming 
the 'Small Magellanic Cloud' extinction curve). We check all three 
sources for archival spectra with broad H/3 an d Ha lines from 
which to estimate the intrinsic reddening E(B-V) (Ward etal.ll987h 
but do not find any distinctive broad components. We therefore pro- 
duce an estimate of the intrinsic reddening by converting from the 
already determined X-ray Ah to a p redicted dust extinct ion E(B- 
V) using a standard dust-to-gas ratio (Maio lino et alj200lh . For fit- 
ting the combined opt ical, UV and X - ray da ta, we identically fol- 
low the approach in Vasu devan et al.l d2009l) . The DISKPN model 
has the following parameters: the inner radius of the disc R ln , the 
maximal temperature in the disc T max and the normalisation TV = 
M BH cos{i) /Dl/3 4 (for black hole mass Mbh, inclination angle 
i, luminosity distance Dl and colour-to-effectiv e temperature ra- 
tio 0) . We assume i = and /3 = 1 (as done in I Vasudevan et all 
2009) and R ln = 6 R g as is appropriate for efficient accretion. 
We fit the optical-UV spectrum with the reddened accretion disc 
model combination to determine T max , which provides an estimate 
of the turnover in the unobserved far-UV. The X-ray (XMM-PN) 
data were fit in conjunction with the optical-UV OM data, retain- 
ing the best-fitting X-ray model above 0.4 keV, giving a full model 
combination of ZDUST(DISKPN)+TBABS*ZTBABS(POWERLAW + 
ZTBABS*(POWERLAW+ ZGAUSS)). After a satisfactory fit was ob- 
tained, all absorptions (intrinsic-to-source and Galactic; both dust 
reddening and neutral gas) were finally set to zero before deter- 
mining the 'primary' (direct) estimate of the bolometric luminosity, 
defined as L bo i = Lo.ooi-aookeV- 

For the second approach, we assemble data from Spitzer 



where available, and supplement it with longer-wavelength Her- 
schel data where available. We use the approach outlined in V10 
and fit a nuclear and ho st galaxy SED template (using the templates 
from lSilva et al.ll2004l) . The nuclear SED templates are generated 
from a combination of mid-IR data and interpolation between dat a 
points using the dust emission model o f lGranato & Danesel ( fl994h . 
and the host galaxy templates are generated using a combination of 
mid-to-far IR data (with th e nuclear contribu tion subtracted), em- 
ploying the GRASIL cod^fl (Silva et al. 1998) for active and quies- 
cent stellar populations to interpolate in the gaps not covered by the 
data. 

In this approach, different regimes of X-ray column density 
Ah are found to produce different IR spectral shapes for the nu- 
clear emission, and different X-ray luminosities are found to pro- 
duce different shapes for the host galaxy emission. Using the X-ray 
determined A^h value, we select the a ppropriate nuclear IR tem- 
plate from the available ones provided in lSilva et al.l ( 120041) , and we 
use the X-ray 2-10 keV luminosity to similarly select the appropri- 
ate host galaxy IR SED template. We fit both templates jointly to 
the available near- and far-IR data. We determine the bolometric 
luminosity by integrating the nuclear IR component between 8 and 
1000/im, multiplying it by appropriate geometry and anisotropy 
corrections, and adding to the 0.5-200 keV absorbed X-ray lumi- 
nosity. 

Th e corrections for geometry and anisotropy are described 
fully in IPozzi et al.l 120071 and V10 and are necessary to recover 
the full reprocessed emission. The geometry correction corrects for 
the absorbing medium only partially covering the central accretion 
disc; it therefore only absorbs and re-emits part of the emission. 
Ba sed on the fraction of absorbed-to-unabsorbed AGN required by 
the lGilli et alj d2007h the covering fraction of the torus is estimated 
to be / ~ 0.67 and therefore the observed nuclear IR emission is 
scaled by fc gC0 m = 1// to recover the primary emission (this value 
for the geo metry scaling ha s been confirmed by the more recent 
analysis of lPozzi et al J|20 loh . 

The anisotropy correction takes into account the fact that the 
nuclear IR spectra of sources with different Ah differ significantly 
in the 1-30/^m band, with higher A^h AGN exhibiting a notable 
deficit in this region due to an isotropic emission from the torus 
(see Fig. 1 of ISilva et alj|2004h . Under the unified model, lower 
Nh sources are more 'face-on' to the observer, so to recover the 
true total emission from higher Nh, more inclined sources, we ap- 
ply a correction to recover the 'missed' flux that escapes from our 
line-of-sight. We adopt the Seyfert 1 SED template as the true spec- 
tral shape of a face-on AGN and normalise all other Nh -group SED 
templates to have the same flux as the Seyfert 1 SED in the 30-100 
pim band (where anisotropy differences are negligible for differ- 
ent Nh groups). The correction factor for each Nh group is then 
computed by dividing the integrated flux of a Seyfert 1 SED tem- 
plate over 1-30 /im (where anisotropy differences are most pro- 
nounced) by the obser ved flux in the sa me band for each of the 
Ah groups (as done bv lPozzi et al.l2007h . The correction factor in- 
creases with Ah and we adopt the same corrections as used in V10, 
i.e. fc all is = 1.3 for 22 < logNH < 24 sources and fc an is = 3.5 
for logNH > 24 sources. We finally add the corrected IR nu- 
clear emission (-kg^ioooum x ^gcom x fcanis) to the emission that 
emerges at hard X-ray wavelengths (i.e., the 'remainder' of the in- 
trinsic emission that gets through the absorption, up to 200 keV) to 
determine the total Lboi- 



2 http://adlibitum.oat.ts.astro.it/silva/ 
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We perform a Monte-Carlo analysis to determine errors at 
each stage, but find that systematic effects dominate and that the 
choice of method to determine Lboi prouces a much more pro- 
nounced variation on estimates than the formal errors would in- 
dicate. We therefore quote ranges for Lboi and derived parameters, 
rather than providing formal errors. All key results from this anal- 
ysis on our three objects are presented in Table [3] 

4.2 NGC 454E 

In the OM images, we identify the eastern source as being respon- 
sib le for the AGN em ission based on line ratio diagnostics using 
the I Johansson! |l988) data and the other considerations discussed 
in i ]3,l | In order to estimate the intrinsic extinction, we use the 
Maiolinoe t"al]d200lh dust-to-gas ratio (Fig. 1 of their paper). Their 
paper shows hints that the ratio log(E(B — V) /Nh) reduces at low 
intrinsic X-ray luminosities, below ~ 10 42 erg/s, with a sharp dis- 
continuity leading to very high values of this ratio for even lower 
luminosities. At the 2-10 keV luminosity of this s ource, there is a 
large u ncertainty in the dust to gas ratio as given bv lMaiolino et alj 
d200ll) . and if we use log(E(B - V) /Nh) = -23.6 to -23.0 we 
find a range of values for the extinction 0.9 < E(B — V) < 3.68. 

This uncertainty does not, however, produce a correspond- 
ingly large uncertainty in the bolometric luminosity. If we use 
E(B - V) = 0.9, fit the combined X-ray and OM data with this 
value fixed into the dust extinction component, and integrate the 
de-reddened SED to determine the bolometric luminosity Lboi, we 
find a value of Lboi = 7.01 x 10 43 erg/s. Using the value of the 
intrinsic X-ray luminosity determined before, we find a bolometric 
correction to X-rays of 41, typical of that for Seyferts. Employing 
the other extremal value of E(B — V) = 3.68 only produces a 
5% change in Lboi and a very similar bolometric correction (38). 
We therefore adopt E(B — V) = 0.9 for all s ubsequent calcula- 
tions. The Eddington ratio estimate, using the I Winter et alj (2009) 
estimate of black hole mass (log Mbh=6.23) is AEdd = 0.32. 

Inspection of Hubble Space Telescope (HST) Wide-Field 
Planetary Camera 2 (WFPC2, filter 814W at 8140A) data for this 
source using the XIMAGE software yields very little hint of any 
strong nuclear point source, as is expected for Seyfert 2 nuclei 
- the emission profile is entirely consistent with diffuse emission 
from the host galaxy. So, we must consider the possibility that the 
OM emission, with much poorer resolution (larger point-spread- 
function, full- width half-maximum of ~ 2") than the HST data 
(FWHM for the WFPC2 ~ 0.1"), represents a strong upper limit 
on the true AGN emission, and may give us a significant overesti- 
mate of the bolometric luminosity. However, the presence of [Ne 
V] high-ionisation emission lines at 14.32 and 24.32/im in the nu- 
clear mid-IR spectra of this object are a clear indication that AGN 
activity is present, and that there must be some emission from 
the c entral engine in th e far-UV at E > 98eV to produce these 
lines dWeaver et al.l201oT) . The severely reddened optical-UV SED 
shape is therefore highly likely to be indicative of dust. 

For our second approach, we fit archival Spit zer spectra 
with t he host-galaxy and nuclear AGN templates from lSilva et alj 
j2004h . We used Spitzer spectroscopy observations in the Short- 
High (SH) and Short-Low (SL) IRS order in the staring mode 
(5.5 < A < 19^im). These spectroscopic apertures have the 
smallest angular resolution in the dispersion direction (3.6-4.7") 
and because of the slit positions in the sky, the SH and SL mod- 
ules only include the nuclear contribution from NGC 454E. On the 
other hand, the extracted fluxes from the Long-High and Long-Low 
IRS order are contaminated with emission from the western com- 



panion galaxy, resulting in a continuum flux approximately twice 
the flux from the smaller slits alone. We therefore restrict our fits to 
the short-wavelength SH and SL spectra. 

We use the intrinsic X-ray luminosity and absorption Nh to 
select the appropriate IR SED templates to fit; for NGC 454E the 
column of ~ 3 x 10 23 cm' 2 implies that the nuclear template 
for log(A r H)=23-24 is most appropriate for fitting to the IR data. 
However, the substantial variation in absorp tion in this sourc e sug- 
gests that the Compton-thick template in ISilva et ID d2004h may 
also be appropriate. However, we find that varying the choice 
of nuclear template used does not produce significantly differ- 
ent results for Lboi. The bolometric luminosity estimate fr om this 
method is 2.4 x 10 43 ergs _1 . The iMarchese et alj d2012l) analy- 
sis of the Suzaku data reveals a significantly different spectrum, 
and if we employ their best-fit model to determine the total 0.5- 
200 keV luminosity for the bolometric luminosity calculation, we 
find Lboi = 2.8 x 10 43 ergs _1 . Gathering together the differ- 
ent estimates of Lboi using the IR SED integration approach and 
using log(MBH) = 6.23, this corresponds to an Eddington ra- 
tio of 0.11 — 0.13, considering all possible X-ray spectral shapes 
(from XMM and Suzaku analyses) and column density estimates. 
The bolometric correction to 2-10 keV using the XMM observa- 
tion is 15, whereas it reduces to 4 if we use the results from the 
Suzaku analysis. Both the 'direct' (optical-to-X-ray) estimation 
and the 'reprocesse d' measures of Lboi yield values less than the 
iFabian et al] d2009h estimate, pushing this object further away from 
the 'forbidden zone' of the Nh — AEdd plot into the 'long-lived 
absorption' region. Notably, the IR-based estimate of Lboi is sig- 
nificantly less than the optical-to-X-ray SED integration estimate, 
indicating that there is much uncertainty in the dereddening correc- 
tion and location of the turnover in the unobservable far-UV part of 
the spectrum. For this source, the optical-to-Xray SED integration 
method may severly over-estimate the bolometric luminosity. 



4.3 2MASXJ03565655-4041453 

For this object, the intrinsic X-ray luminosity of 3.5 x 
10 4 3 erg/s implies E(B - V)/Nn is ~ 10~ 23 magcm 2 from 
the iMaiolino et alj d200lh plot, yielding an intrinsic extinction of 
E( B — V) = 0-71. F itting the broad-band SED as before, using 
the lWinter et alj J2009) estimate of black hole mass to constrain the 
normalisation of the DISKPN model, and hence, the position of the 
unobservable 'turnover' in the big blue bump, we find a bolometric 
luminosity of 1.48 x 10 45 erg/s. For this object, we also have an es- 
timate of the black hole mass using the V09/V10 refined method of 
dealing with the extended and nuclear K-band emission, and using 
this mass instead in the fit yields Lboi = 1-83 x 10 45 erg/s, due to 
the different predicted turnover in the spectrum obtained using this 
mass estimate. The Eddington ratios obtained using the two mass 
estimates are 0.03 (W09 mass) and 0.06 (V10 mass), with bolo- 
metric corrections of 42 (W09 mass, log M B h = 8.64) or 52 (V10 
mass, log Mbh = 8.38). Combining results for the Eddington ra- 
tio and column density, we again find the object pushed away from 
the 'forbidden region' into the 'long-lived absorption' regime. We 
cannot find any usable archival IR data for this source to estimate 
Lboi using the reprocessed emission. 



4.4 XSSJ05054-2348 

The intrinsic X-ray luminosity of ~ 10 43 erg/s imp l ies E (B — 
V)/Nh is ~ 10~ mag cm from the lMaiolino et alj d200ll) plot, 
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yielding an intrinsic extinction of E(B — V) = 0.1. Fi tting 
the broad-band SED as before, using the I Winter etail (2009) es- 
timate of black hole mass, we find a bolometric luminosity of 
2.59 x 10 44 erg/s. We have another estimate of the black hole 
mass from V10 and using this mass instead in the fit yields Lboi = 
2.55 x 10 44 erg/s. The Eddington ratios obtained using the two 
mass estimates are 0.06 (W09 mass) and 0.04 (V10 mass), with 
bolometric corrections of 8.6 (W09 mass, logMeH = 7.53) or 
8.4 (V10 mass, log M BH = 7.67). 

We also produce an IR-based estimate of the bolometric lumi- 
nosity using archival Spitzer and Herschel data. There is no con- 
fusion between the source of interest and any background sources 
in either dataset. We used Spitzer spectroscopy observations in the 
Short-Low and Long-Low IRS order in the staring mode (5.5 < 
A < 36/im). We follow the same procedure to determine Lboi as 
for NGC 454E, described in S |4.2| fitting host and nuclear compo- 
nents, and determining Lboi by adding the absorbed 0.5-200 keV 
luminosity to the nuclear IR luminosity from 8-1000/im with ap- 
propriate corrections for geometry and anisotropy (see ^4. 1 1 >. The 
bolometric luminosity from this method is 4.4 x 10 44 erg s _1 . This 
is a factor of 1.7 larger than the estimate from de-reddening the 
optical-UV SED and integrating, emphasising the uncertainty in 
the de-reddening of the optical-UV spectrum. The corresponding 
2-10 keV bolometric correction is 14.4 from this method, and the 
Eddington ratio estimates range between 0.07-0.12 (the range is 
from using the V10 mass or the W09 mass estimates). Combining 
results for the Eddington ratio and column density, we again find 
the object pushed away from the 'forbidden region' into the iong- 
lived absorp tion' regime, compa red to the position of this object as 
presented in Fab ian et"al] d2009l) . 



5 CONSTRAINTS ON Lbol AND A E dd FROM MID-IR 
HIGH-IONISATION LINES 

The work of Melende z et alj d201 ll) provides a constraint on the 
broad-band SED shape and bolometric luminosity. They generated 
photoionisation models to characterize the small dispersion found 
in correlations between high-ionisation mid-infrared emission-line 
ratios in AGN in order to predict the slope connecting the intrin- 
sic emission at 1 Ry to the intrinsic emission at 1 keV, providing 
a very useful constraint on the unobservable extreme-UV (EUV) 
emission. Other works have investigated the relationship between 
other high-ionisation IR lin es and the unobserva ble UV emission 
in individual objec t s (e.g ., lMullanev et al.ll2009l for Ark 564), but 
the lMelendez et al] J20 1 lh study constitutes the first attempt to con- 
strain the unobserved EUV continuum in the statistically significant 
BAT AGN sample, and provides a useful empirical tool for connect- 
ing the readily observable 1 keV emission to the unobservable EUV 
emission. This allows us to estimate whether our bolometric lumi- 
nosity estimates in the previous sections are sensible, using simple 
physical constraints on the UV energy required to produce mid IR 
emission lines. For NGC 454E, we are able to obtain the mid-IR 
line ratios from IWeaver et all d2010t) and predict an energy index 
a — 1.7 for the power-law joining the 13.6 eV and 1 keV emission 
(for F v oc v~ a ). If we extrapolate from the intrinsic 1 keV emis- 
sion using this spectral index, we see from Fig.|7]that the emission 
connects well with the predicted 13.6 eV emission from our accre- 
tion disc model fit using the de-reddened DISKPN model. Whilst the 
IR-based estimate of Lboi is significantly lower than the estimate 
from the UV DISKPN fit, the high-ionisation line ratios would im- 



ply that the Lboi estimate from the UV model fit is consistent with 
the physical constraints from the high-ionisation mid-IR lines. 

We do not have high-ionisation lines available in the mid-IR 
spe ctrum for XSS J05054 -2348 or 2MASX 103565655-4041453, 
but lMelendez et all d201ll) identify that the range of 13.6 eV-1 keV 
energy indices is 1.5-2.0 for the BAT AGN sample. For XSS 
J05054-2348, if we use the minimal index 1.5 which produces the 
minimal possible emission at 13.6 eV, we find that it far exceeds 
that predicted from the DISKPN big blue bump fit (Fig. [9}. In this 
source, optical-to-X-ray SED fitting is likely to severely underes- 
timate Lboi, and the IR-based estimate is more likely correct. For 
2MASX J03565655-4041453, we show the range of possible val- 
ues for the 13.6 eV emission using 1.5 < a < 2.0, and find that the 
de-reddened DISKPN model fit lies well within these limits (Fig.[8]l. 



6 DISCUSSION 

6.1 The Effective Eddington Limit Paradigm 

The refined estimates of column density and Eddington ratio using 
the high-quality XMM-Newton data suggest that all three objects 
under consideration here move away from the 'forbidden zone' 
where short-lived absorption and outflows are expected. These re- 
sults are collated and summarised in Fig.[l0] along with the results 
for the re st of the 9-month BA T catalogue that were previously pre- 
sented in Fabia n et al. ||200i). We again draw the reader's attention 
to the other source (apart from the three studied in this paper) which 
was already known to lie well into the forbidden region, MCG- 
05-23-16, which is known to hav e a complex and variable warm 
absorber (e.g. iReeves et alj |2007[) and a possible variable, high- 
ionisation, high-velocity outflow dBraito et al.l2007t) . Other objects 
near the bo undary are NGC 526 A, which is found to have a warm 
absorber by iNandra et all d2007l); NGC 3783, which has a known 
warm absorber from iKaspi et all d200lh : and IC 4329 A, where a 
prominent dust lane visible in opti cal images may be resp onsible 
for a large part of the A^h seen, and Markow itz et alj d2006l) find an 
outflow . Taken in conjunction with the discovery of IWinter et alj 
d2012h (Fig. 13 of their paper) that warm absorber strength is high- 
est near the Effective Eddington Limit, the evidence increasingly 
points to a scenario where the Effective Eddington Limit is of key 
importance in understanding the presence and nature of dusty gas 
absorption in AGN. 

Returning to the objects newly-observed with XMM studied 
in this paper, we note that in the case of NGC 454E, significant 
absorption variability betw een the XRT, XM M and Suzaku obser- 
vations is observed, and the Marchese et alj d2012l) analysis of this 
XMM dataset comments on a distinct excess between 5-7 keV that 
could be indicative of an ionised absorber with significant column 
(6 x 10 23 cm -2 ), an ionisation parameter of log(£)=3.6 with a tur- 
bulent velocity of 300 kms -1 , emerging from within the broad- 
line region. However, we notice that the object has moved from the 
'forbidden region' (in the Swift/XKF observation) to the 'long-lived 
absorption' region with a Compton-thick column density (in the 
Suzaku observation) to a less-absorbed state also in the 'long-lived 
absorption' region, with an ionised absorber (in the XMM observa- 
tion). It is difficult to connect the properties of this ionised absorber 
directly with the Effective Eddington limit paradigm because of this 
unusual trajectory in parameter space. 

The class of AGN in which absorption variability is present 
('changing l ook AGN') have been discussed extensively in the lit - 
erature (e.g.. lRisaliti et al.ll2010l [20091 120021: [Bianchi et alj|2005h . 
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AGN 


r(l) 
L ha\ 

(1) 




r(2) 
^bol 

(2) 


ft '2-10keV 

(3) 


(2) 

K 2-10keV 

(4) 


A (1) 
A Edd 

(5) 




x (2) 
A Edd 

(6) 


NGC 454E 


7.01 




2.40 


41 


15 


0.32 




0.11 


2MASX J03565655-4041453 


148 - 


183 




42 - 52 




0.03 - 


0.06 




XSS J05054-2348 


25.5 - 


- 25.9 


44.0 


8.4-8.6 


14.4 


0.04 - 


0.06 


0.07 - 0.12 



Table 3. SED fit results using XMM data in X-rays. All luminosities are in units of 10 4s erg s — 1 , (1): bolometric luminosity from the optical-to-X-ray SED 
integration method (method 1). (2): bolometric luminosity from summing reprocessed nuclear IR emission and hard X-ray emission (method 2). (3): X-ray 
bolometric correction, -Lbol/^2— lOkeV using method 1 (4): X-ray bolometric correction from method 2. (5): Eddington ratio, ^bol/^Edd. from method 
1. (6): Eddington ratio from method (2). For objects where two estimates of the black hole mass (W09 and V10) are available, we quote ranges in derived 
quantities arising from using two different black hole masses. 




Figure 10. Changes in absorbing column and Eddington ratio for the three objects with new XMM data analysed in this paper. Small circles represent objects 
with fits to archival ASCA or XMM-Newton data from Winter et al. (2009), and squares represent objects for which only XRT data were available at the 
time that W09 was published. Black filled circles or squares depict objects for which the best fit was a 'complex' spectrum, and empt y circles or squares 
represent those with 'simple' model fits. Red circles show the positions in the Nn — ^Edd of the three objects with new data as given in iFabian etalJ l 2009) 
using XRT data only; lines connect these circles to the new positions of these objects as calculated from the new XMM data. Blue stars show the extremes 
of the Eddington ratios that may be obtained using different black hole mass estimates, intrinsic reddening corrections or Lbol determination methods (i.e. 
optical-to-X-ray SED integration vs. using the IR reprocessed emission). In all cases, the full range of Eddington ratios obtained, if combined with the new 
N-r estimates, place all three objects within the 'long-lived absorption' region. 



These AGN could potentially represent cases where radiation pres- 
sure acts to change the absorption on longer timescales, without 
necessarily manifesting as outflows. 



7 SUMMARY AND CONCLUSIONS 

Using high-quality XMM Newton spectra, we present new spectra 
for three objects from the 9-month Swift/BAT catalogue that were 
identified as being close to the Effective Eddington Limit for dusty 
gas. The XMM spectra reveal that the best-fit model is a double- 
power law (partial covering absorption) for all three objects. In each 
case, the scattered component in the second (low-energy) power- 



law for all three objects is ~1 per cent of the primary power-law 
(emergent at high energies in the BAT spectra), consistent with sim- 
ple electron scattering, not requiring any particularly sophisticated 
geometry for the absorption. 

These spectra from XMM uncover hidden complexity in the 
spectrum for two of our objects, 2MASXJ03565655-4041453 and 
XSSJ05054-2348, for which the previously obtaine d XRT data 
were adequately fit by simple power-laws in the Win ter et al.l 
(2009) analysis. This highlights the usefulness of XMM in un- 
earthing nuances in the spectra of objects, providing more infor- 
mation about the environment in these objects. 

Taking into account uncertainties in the black hole mass esti- 
mates and estimates of intrinsic reddening, we have produced es- 



© 0000 RAS, MNRAS 000, 000-000 



12 R.V. Vasudevan et al. 



timates of the bolometric luminosities and Eddington ratios using 
the new XMM data. We employ the XMM OM data to generate 
broad-band SEDs for the three objects and attempt to recover their 
bolometric luminosities by applying a first-order reddening cor- 
rection using a dust-to-gas conve rsion factor from the X-ray col- 
umn density (Maiolino et al. 200 1]). We also use the re-processed 
infrared emission, as traced by Spitzer and Herschel observations 
where available, to obtain independent estimates of the bolometric 
luminosity and compare them with those obtained from optical- 
to-X-ray SEP integration. We us e estimates of the black hole mass 
using the lVasudevan et ai1 ([2~009) method if available; in the case of 
NGC 454E we only have one M BH estimate (from W09). The Ed- 
dington ratios determined all lie further into the 'long-lived absorp- 
tion' regime than found previously. NGC 454E exhibits the highest 
variability in the absorption. This AGN also exhibits the highest 
Eddington ratio of the three; the variable absorption places it in the 
'changing-look AGN' category (Risaliti et al. 20020, so it may be 
the case that the changing absorption is driven by radiation pressure 
from the central engine. However, there is little evidence so far to 
suggest that changing-look AGN as a class are consistently located 
at the Effective Eddington limit. 

This study reinforces the need to use high-quality spec- 
tra to determine the nature and amount of absorption in AGN, 
and re-emphasises the need for good mass estimates to deter- 
mine Eddington ratios. The recent comprehensive study 58-month 
BAT catalogue AGN in the Northern Galactic Cap (b > 50°, 
IVasudevan et alj |2013) presents refined, systematically determined 
JVh estimates for the latest, deepest version of the BAT AGN cata- 
logue. If coupled with carefully deter mined black hole ma sses, this 
could provide a useful revision of the iFabian et al.1 j2009t) work on 
the earlier BAT catalogue, offering new insights on the interplay 
between radiation pressure and absorption in local AGN. 
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